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HYDROGEN  SULFIDE  EFFECT  ON  HYDROGEN  ENTRY  INTO  IRON 

-  A  MECHANISTIC  STUDY 


Rajan  N.  Iyer,  Izumi  Takeuchi*,  Mchrooz  Zamanzadeh**  and 

Howard  W.  Pickering 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

ABSTRACT 

The  recently  developed  I-P-Z  model  is  modified  in  order  to  analyze  the  observed  enhanced 
permeation  of  hydrogen  that  occurs  in  the  presence  of  hydrogen  sulfide  during  cathodic  hydrogen 
charging  of  iron.  The  modification  accounts  for  the  fact  that  the  energy  of  adsorption  becomes 
coverage  dependent  at  the  higher  coverages  and  affects  the  hydrogen  evolution  reaction  (h.e.r.)  in  the 
presence  of  H2S.  Charging  experiments  were  performed  on  Ferrovac  E-Iron  membranes  0.5  mm 
thkk  using  a  Dcvanathan-Stachurski  cell  in  deaerated,  pre-dectrolyzed  solutions  made  from  O.IM 
HCIO4  and  0.1M  NaClC>4  with  pH  values  of  1  and  2.  The  transfer  coefficient, ....  exchange  current 
density,  i0,  thickness-dependent  absorption-adsorption  rate  constant,  k”,  recombination  rate  constant, 
kj.surface  hydrogen  coverage,  8h.  and  discharge  rate  constant,  k°,  were  obtained  by  application  of 
the  model  to  the  experimental  results.  As  a  result,  the  role  of  HoS  has  been  clarified.  While  0H  is 
increased  in  the  presence  of  H9S,  the  ovcrpotential,  rj»  is  decreased  consistent  with  an  observed 
increase  in  a,  and  the  increased  H  entry  is  found  to  be  the  result  of  a  decreased  k3  as  well  as  the 
increased  a.  In  addition,  a  very  important  relationship  has  been  derived  that  v>ill  enable  the  calculation 
of  the  absorption,  kabS,  and  adsorption,  kiKjs,  rate  constants  from  the  electrochemical  permeation  results 
for  different  membrane  thicknesses. 

*  Present  address:  Kashima  Steel  Works,  Sumitomo  Metals  Industries,  Ltd.,  Kashima,  Japan. 

**  Present  address:  Professional  Services  Industries  -  P  T  L  Division,  850  Poplar  Street,  Pittsburgh, 
PA  15220. 


INTRODUCTION 


Since  the  failure  of  steels  in  wet  hydrogen  sulfide  (H2S)  is  well  known  to  occur  in  many 

technological  situations,  especially  in  oil  and  gas  field  applications  many  studies  have  been  carried 
out  in  order  to  better  define  and  understand  this  problem.  It  is  well  known  that  in  the  presence  of 
H2S,  the  concentration  of  hj  d^ogen  in  steels  is  increased  sharply  U-6).  ft  js  ajso  ll0W  generally 

accepted  that  when  the  concentration  of  hydrogen  in  a  metal  reaches  a  critical  level,  cracking  can  occur 
(7). 

It  has  often  been  suggested  based  on  a  measured  increase  in  hydrogen  uptake,  that  I-I2S  has 

a  "poisoning"  effect  on  the  hydrogen  evolution  reaction  (n.e.r.).  Contrary  arguments  based  on 
polarization  and  permeation  experiments  have  provided  different  mechanisms  13,4,5,8)  for  (j,c  roic  0f 
H2S  in  the  entry  of  hydrogen  into  iron  and  steel.  There  is  also  a  mechanistic  dilemma  on  the 
significant  effect  that  H2S  has  on  the  kinetics  of  the  h.e.r.  and  hence  on  hydrogen  absorption. 

There  arc  at  least  four  main  proposals  of  the  mechanism  of  the  enhancement  of  hydrogen 
absorption  by  iron  due  to  PI2S.  The  first  of  these,  which  was  developed  by  Iofa  and  Katn  ® 


considers  the  following  sequence  of  reactions: 

Step  1:  MS'  =  (HS-)ads 

Step  2:  <HS-),,ds  +  H30+  =  (H-S-H)ads  +  H20 

Step  3:  (H-S-H)ads  +  e'  — >  (HS')ads  +  Hads 


In  this  model,  steps  (1)  and  (2)  are  quasi  reversible  and  step  (3)  is  considered  to  be  rate  limiting.  On 
the  basis  of  their  experimental  data  in  0.1N  H^O^  Iofa  and  Kam  concluded  that  the  HS'  anion 

accelerates  the  discharge  of  H+  ion.  However,  their  interpretations  were  incomplete  in  the  sc  ;  that 
the  reaction  parameters  remained  unknown  and  in  fact,  there  is  a  far  higher  concentration  of  H?S 

compared  to  HS*  in  acid  solutions  as  shown  in  Figure  1.  So,  Kawashima  ct  al  proposed  the 


1  1 


ja* 


following  reaction  sequence: 

Step  1:  H2sads  +  e'  *  H2sads~ 

Step  2:  H2Sa(]s“  +  Hads+  H2S . Hads 

Step  3:  ^2^ . ^ads  ^  ^2^ads  +  ^ads 
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They  showed  that  IbS  did  not  change  the  Tafel  slope,  but  did  decrease  the  cathodic  overpotential  in 
acid  solutions.  In  their  model,  H2S  works  as  a  bridge  to  transfer  the  electron  from  the  iron  surface  to 
the  hydrogen  ion. 

The  poisoning  model  of  Bockris  et  al  attributes  the  increase  in  hydrogen  absorption  to  a 
lowering  of  the  M-Ha(js  bond  energy,  w.  Lowering  of  w  increases  the  activation  energy  for  the 

discharge  process  of  the  h.e.r.  Thus,  in  systems  where  this  model  is  applicable,  the  hydrogen 

overvoltage  would  increase  with  increasing  concentration  of  the  poison.  However,  it  is  now  well 
established  that  H2S  decreases  the  hydrogen  overpotential  (3, 4, 5, 8)  anc^  thUSj  ,hc  poisoning  model  is 

not  applicable  in  this  case. 

Bolmer  ^  ^  proposed  an  overall  reaction  model: 

2  I-bS  +  2  e*  -»  H2  +  2  HS- 

and  rationalized  that  at  currents  below  the  diffusion  limited  current,  the  hydrogen  evolution  mechanism 
is  a  function  of  the  HS"  concentrauon  and/or  the  H2S/HS‘  ntuo.  However,  his 

analysis  was  incomplete  in  that  the  reaction  parameters  for  the  h.e.r.  and  hydrogen  absorption 
remained  unknown. 

In  the  studies  cited  above,  the  experimental  data  were  not  sufficient  and  the  interpretations  were 
based  mostly  on  experiments  which  were  not  well  defined.  Also,  up  to  the  present  time  and  because 
of  its  complexity  there  seems  to  be  no  general  model  which  can  rationalize  the  significant  effect  of  H?S 

on  the  h.e.r. 

This  paper  reports  the  results  of  experimental  studies  directed  at  evaluating  the  effect  of  H^S  on 

the  hydrogen  permeation  and  evolution  reactions.  It  analyzes  these  results  with  a  newly  developed 
general  model  of  hydrogen  entry  into  metals  ^2,13)  jn  orcjer  t0  gajn  ,m  jnsjg|u  jni0  L|ie  permeation 
mechanism  and  to  arrive  at  otherwise  difficult-to-obtain  parameters  such  as  the  surface  coverage  and 
the  individual  rate  constants. 
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EXPERIMENTAL 

A  Dcvanathan-Stachurski  cell 1  ^  was  used  to  measure  die  permeability  of  hydrogen  through 
iron  membranes  as  a  function  of  the  charging  current,  pH  and  IhS  concentration  in  the  charging 

solution.  The  apparatus  shown  in  Figure  2  consists  of  two  identical  electrolytic  cells  separated  by  the 
iron  membrane.  One  of  the  membrane  surfaces  (input  surface)  was  galvanostatically  charged  with 
hydrogen,  using  a  Pi  counter  electrode.  The  potential  of  the  input  surface  was  measured  with  a 
saturated  calomel  reference  electrode.  The  other  membrane  surface  (exit  surface)  was  coated  with  Pd 
and  held  in  0.1M  NaOh  at  a  potential  of  150  mV  (versus  a  HgO/Hg  reference  electrode)  which 
oxidized  all  of  the  hydrogen  diffusing  through  .he  membrane;  this  oxidation  current  provided  a 
measure  of  the  hydrogen  permeation  flux. 

The  membrane  was  Fcrrovac  E-Iron,  annealed  m  evacuated  capsules  at  930°C  for  two  hours  and 
furnace  cooled.  The  final  preparation  involved  polishing  with  600  emery  paper  and  degreasing  with  a 
mixture  of  benzene  and  acetone.  Tire  final  thickness  of  the  iron  specimen  was 
0.5mm  and  this  iron  membrane  was  fixed  between  the  two  cells,  tightly  sealing  them  with  rubber  O- 

rings.  Details  of  the  cell  design,  circuitry,  edge  effects,  and  special  procedures  can  be  seen  elsewhere 
(15). 

Mixtures  of  0.  i  HCIO4  and  0.1M  NaC104  were  used  to  make  up  the  charging  solutions  with 

pH  values  of  1  and  2.  The  solutions  were  prc-electrolyzed  for  72  hours  under  an  atmosphere  of 
prepurified  nitrogen.  To  make  up  solutions  containing  H2S,  Irst  a  master  solution  0IT-I2S  was 
prepared  by  passing  H2S  gas  of  99.67o  purity  through  the  deaerated,  pre-clectrolyzed  0.1M  HCIO4 
solution.  The  concentration  of  H2S  in  this  master  solution  was  measured  with  an  ion  selective 
electrode  and  found  to  reach  5.2x10’“  M,  after  two  hours  of  H^S  flow  through  the  solution.  The 
master  solution  of  H2S  was  mixed  with  the  charging  solution  in  order  to  give  H2S  concentrations  of 
10*5  to  10’3  M,  while  maintaining  the  ionic  strength  constant. 

In  order  to  obtain  reproducibility  of  the  permeation  data,  several  different  precharging  procedures 
were  evaluated.  The  procedure  that  gave  good  reproducibility  consisted  of  the  following  steps.  The 
charging  solution  was  pre-electrolyzed  in  a  pre-electrolysis  cell  with  Pt  electrodes  at  2  mA  for  more 
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limn  72  hours  prior  10  starting  the  experiment.  Then,  the  iron  membrane  was  precharged  ai  l  mA  for 
50  minutes  and  then  ai  0.5  mA  for  3  hours  {during  this  time,  the  charging  potential  moved  in  the 
negative  direction  as  the  permeation  current  increased).  Following  this  step,  the  cell  \\  as  drained  and 
filled  again  with  fresh  pre-eleetroiyzed  solution.  Subsequently,  the  polarisation  and  permeation  data 

were  recorded  proceeding  from  a  starting  higher  current  (=  1mA  cm*")  to  lower  eunrems.  To  perform 
experiments  with  H2S  ,  the  charging  solution  was  changed  to  the  H2S  pre-mixed  solution  and  the 

same  sequence  was  repeated 

RESULTS 

The  cathodic  polarization  curves  for  different  H2S  concentrations  (c„iS)  in  the  charging  solution 

at  two  different  (acidic)  pH  values  are  given  in  Figures  3  (a)  and  (b).  It  is  easily  observed  from  these 
plots  that  the  hydrogen  overvoltage  (T])  decreases  asymptotically  with  increasing  H2S  concentration. 

Figures  4  (a)  and  (b)  display  the  steady  state  permeation  current  (i„)  as  a  function  of  the  hydrogen 
overvoltage  (T|),  for  the  pHl  and  pH2  charging  solutions.  From  these  figures,  it  is  obvious  that  i,., 
increases  substantially  in  the  presence  of  H2S.  They  also  show  that  while  all  of  the  In  ic  vs  T|  and  In 
L  vs  r|  plots  in  the  non  H2S-coninining  solutions  are  linear,  not  all  are  linear  in  the  presence  of  H2S. 

In  order  to  understand  the  mechanism  and  quantitatively  analyze  the  partitioning  between  the 
hydrogen  absorption  and  hydrogen  evolution  reactions,  the  recently  developed  I-P-Z  model  U3)  js 
applied  here.  According  to  this  model,  the  i„  vs  V ir  plot  (where  ir  =  ic-i,„  is  the  hydrogen 
recombinauon  current  density)  and  fiiq  (=iccaarl)  vs  i*,  plot  (where  a  is  the  h.e.r.  transfer  coefficient 
and  a  =  F/RT  =  38.94  V**)  will  be  linear  if  the  mechanism  of  the  h.e.r.  is  a  non-activated  (i.e.,  the 
energy  of  adsorption  is  not  coverage  dependent)  discharge-chemical  recombination  process.  Figure 
5(a)  shows  that  the  reladonship  between  iM  and  Vir  is  linear  at  pH=l,  for  solutions  with  and  without 
H2S.  Figure  5(b)  shows  that  the  iM  vs  Vir  plot  is  linear  at  pH=2  in  the  absence  of  H2S  and  also  when 
0^5=2.5x10°  M  whereas  it  is  nonlinear  for  the  higher  H2S  concentrations.  The  corresponding  flTj 
vs  i„  plots  are  given  in  Figures  6  (a)  &  (b).  Although  there  is  a  fair  amount  of  scatter  in  Figure  6(a), 
mainly  due  to  the  lack  of  knowledge  of  the  exact  a  values,  the  fj^  vs.  i„  plot  is  found  to  be  linear 
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for  the  lower  H2S  containing  pH  1  solutions.  On  the  other  hand,  for  the  pi-12  solutions  (Figure  6(b)), 
the  vs.  i„  plot  is  found  to  be  linear  for  Cj1jiS=0,  but  nonlinear  for  Cj(jiS=2.5:<1Q*^  M.  Thus,  for  the 
higher  pH  and  higher  I  hS  containing  solutions,  deviations  occur  in  the  square-root  relationship 
(Figure  5(b))  or  the  fj^-i^,  plot  becomes  nonlinear  (Figure  6(b)),  or  both. 

ANALYSIS  AND  DISCUSSION 

1 ,  Dam  Analysis  Using  the  Basic  I-P-Z  Model 

Since  the  i„  vs  Vi,,  and  fj^  vs  i„  plots  are  linear  for  the  pHl  solutions  except  at  the  highest  HjS 
concentration,  including  cjUS=0  (Figures  5  (a)  and  6  (a) )  as  well  as  for  the  pH2  solution  at  c,{iS  =  0 
(Figures  5  (b)  and  6  (b)),  the  basic  I-P-Z  model  U3)  can  be  applied  for  these  cases.  The  following 
equations  represent  the  fluxes  or  currents  of  the  discharge,  recombination  and  permeation  reactions 


(symbols  are  defined  at  the  end  of  the  text)  U3); 

ic  =  Fki(l-eH)c-ac°l  =  i0'  tt-0H)«'a“n  0) 

where  i0'  =  F  kj°  3j^+c‘aa^  ^ 

ir  =  Fk3eH2  (2) 

ioo  =  F  (D]/L)  cs  =  cs/b  (3) 

®H=  cs/k"  (3a) 

^  =  ^abs^'ads  +  (Uj/L))  (ab) 

And  from  these  equations,  the  following  relationships  have  been  obtained! 13): 

i„«  r/b)  (F  fc3)-°-5  Vir  (4) 

ic  eaaTl  =  -  (b  i0'/k")  i„  +  i0'  (5) 


Thus,  from  Figures  5  and  6  (in  essence,  by  regression  analysis  of  Equations  (4)  and  (5))  and  Eq. 
(1),  all  of  the  following  constants  are  computed:  the  exchange  current  density,  i0  (=i0'  (l-9e));  the 
discharge  reaction  rate  constant,  k°;  the  recombination  reaction  rate  constant,  and  the  thickness 
dependent  absorption-adsorption  constant,  k".  For  details  of  the  computation  procedure,  see 
elsewhere  U^). 
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2.  Modified  l-P-Z  Model  Incorporating  Frtunkin-Temkin  Correction  for  Activated  Reactions 

For  pH  =  2,  [he  presence  of  H2S  causes  a  deviation  from  a  linear  i„  *<  \!ir  behavior  t  Figure 
5(b))  or  linear  f^ «  i„  behavior  (Figure  6(b)).  One  possibility  is  the  adsorption  of  HS*  but  this 

would  require  that  the  pH  in  the  interface  region  rise  above  -4.5  since  it  is  clear  from  Figure  1  that  the 
dissociation  of  H2S  is  negligible  in  solutions  of  pH  £4.5.  Another  is  that  H2S*  can  be  formed  CI6,17) 

by  the  following  reaction: 

H,S+e-  — »H2$ads  (6) 

having  a  high  reaction  rate  constant.  This  fast  side  reaction  occurring  at  the  cathodic  metal  surface  can 
have  three  different  effects: 

(1)  H2S  can  function  as  a  bridge  for  hydrogen  discharge  through  the  H2S*  species,  similar  to  the 
proposal  by  Kawashima  et  al  $); 

HoSads  +  H30+  ->  (H....H2S)  +  H20  (7) 

(H....H2S)  +  H20  ->  (H....H2S)ads  +  H20  Hads  +  H2S  +  H20  (8) 

where  H2S”  ^  discharges  the  proton  (Eq.  (7)}  and  becomes  the  stable  H2S  molecule  (Eq.  (8)).  This 

step  can  decrease  the  overvoltage  for  the  discharge  reaction,  in  accord  with  the  experimental 
observations  (Figures  3  (a)  and  (b)). 

(2)  Besides  the  above  effect,  the  recombination  reaction  can  also  be  restricted  effectively  by  the 
(H...H2S)  intermediate  of  this  side  reaction  slowing  the  diffusion  of  hydrogen  adatoms  and/or 
blocking  the  hydrogen  recombination  sites  on  the  metallic  surface. 

(3)  At  higher  pH,  in  view  of  the  higher  [H2S]/[H+}.ratio,  chemical  desorption  is  less  strongly  favored; 
i.e.,  effect  2  could  lead  to  a  change  in  the  mechanism  from  a  primarily  chemical  desorption  to  an 
electrochemical  desorption  mechanism. 

The  adsorption  of  HS*  on  me  surface  can  produce  the  same  result  as  the  (H....H2S)  species  in 
effects  (2)  and  (3)  but  HS*  cannot  match  H2S*  in  facilitating  the  proton  discharge  step  (effect  (1)) 
because  HS,  unlike  H2S,  is  not  stable  so  the  bridge  mechanism  is  the  more  likely  one.  The  above 
effects  (2)  and  (3)  indicate  activation  of  the  hydrogen  evolution  reaction  (h.e.r.)  and  under  these 
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conditions,  the  Frumkin-Temkin  (F-T)  corrections  have  to  be  applied  to  the  discharge  and 
recombination  currents  ^);  (hus,  Equations  CD  and  t2)  will  be  modified  as  follows: 

ic  =  V  U-0HU':ianfofeH  (9) 

ir  =  Fk30H2c2afeH  (10) 

where  f = y/RX  and  y  is  the  gradient  of  the  apparent  standard  free  energy  of  adsorption  with  coverage. 
Although  not  considered  in  Eqs.  (8)  and  (9),  the  coverage  by  HjS*  may  also  be  important  as  far  as  site 
availability  and  activation  of  the  reactions  are  concerned.  In  order  to  determine  the  effect  of  H2S,  the 
relationship  between  8h  and  Gh^S*  needs  to  be  quantified  by  considering  the  change  in  0H  by 

additions  of  H2S  to  the  solution.  The  implicit  assumption  in  the  present  model  is  that  this  HoS  effect 
can  be  neglected  as  if  the  H2S*,  in  acting  as  a  bridge  to  discharge  the  proton,  will  not  occupy  a  surface 
site  for  a  prolonged  duration.  A  typical  value  of  f  is  4  to  5  for  hydrogen  coverages  ^°);  f  =  4.5  will  be 


used  here. 

From  Equations  (3)  and  (3a)  we  obtain, 

8H  =  biJk"  (11) 

Inserting  Equation  (1 1)  into  Equation  (10),  it  can  be  shown  that 

In  Wij/iJ  =  (afb/k")L  +  In  {b(Fk3)°-5/k"}  (12) 

Substituting  0pj  from  Equation  (1 1)  into  Equation  (9),  one  obtains 
In  (ic  c^^/d-bi  Jk")]  «*  -  a  a  n  +  In  (i0'  )  (13) 

or 

ln(fic,iJ  =  *aall  +  In(i0')  04) 

where 

f;c.i-.  =  ‘e  e<c‘fbi~/k'  Vo-bi^k")  (14a) 


Thus,  Equation  (12)  tells  us  that  will  not  be  linearly  related  to  Vi,.,  when  the  h.e.r.  is  activated  (f  > 
0),  but  that  i„  will  be  linearly  related  to  In  0'/ir/ieo).  Also,  according  to  Equation  (14),  In  (f,  j^)  will 

be  linearly  related  to  r\  when  the  h.e.r.  is  activated. 
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3.  Data  Analysis  by  the  Application  of  F-T  Corrected  l-P-Z  Model 

Figure  7  shows  that  the  plots  of  In  (Vij/i.J  vs  i„,  at  pH  =  2  sire  linear  in  accord  with  Equation 
( 12)  for  c^alO’ M  and  c,JoS=3xlO'^  M.  But  for  the  lowest  HoS  concentration,  c,ilS=2.5xlO*^  M, 

the  plot  may  not  be  completely  linear  and  is  consistent  with  the  h.e.r.  being  only  partially  activated,  as 
is  already  indicated  above  in  the  i„  vs  Vir  plot  (Figure  5(b))  which  is  linear  and  in  the  f^  vs  i„  plot 
(Figure  6(b))  which  is  nonlinear.  Extended  analysis  showed  that  the  data  for  cHoS=2.5xlO'->  M  are 

belter  analyzed  with  the  F-T  corrected  I-P-Z  model. 

The  value  of  a  is  determined  from  Equations  (12)  and  (14)  by  an  iterative  procedure.  Using  an 
initial  value  of  a  =  0.5  and  the  k"  obtained  from  the  slope  of  the  C\/ir  /i  J  vs.  iM  plot  (Eq.  (12)),  In 
tfj  :)  is  known  for  various  values  of  ic  and  iM  and  as  a  function  of  q.  Thus,  a  regression  analysis 
of  Eq.  (14)  can  be  done  to  obtain  a  new  value  of  a.  Using  this  new  value  of  a,  a  new  value  of  k"  is 
obtained  (from  Eq.  (12))  and  a  regression  analysis  is  again  performed  on  Eq.  (14).  This  procedure  is 
repeated  until  the  a  obtained  from  the  slope  of  the  In  (fj  |^)  vs  q  plot  (Eq.  (14))  converges  to  a  fixed 
value.  Then,  the  In  (fj„tj^)  vs  q  functions  arc  plotted  (Figure  8),  and,  besides  a,  k3,  k"  and  iQ'  are 

determined  from  the  regression  analysis  of  Equations  (12)  and  (14). 

Having  determined  k"  for  all  combinations  of  pH  and  Cj^s  (ic.  with  and  without  the  F-T 

correction  to  the  I-P-Z  model),  0pj  (calculated  using  Equation  (11))  vs  q  can  be  plotted  and  these  are 
shown  in  Figures  9  (a)  and  (b).  These  plots  demonstrate  that  HjS  increases  the  hydrogen  surface 
coverage  (0^)  quite  significantly  and  that  the  increase  is  asymptotic  with  increasing  c^.  In  fact,  at 

pH  =  2  (Figure  9  (b)),  the  0j^  values  increase  significantly  with  small  additions  of  H->S,  but  do  not 
increase  any  further  with  higher  cHoS  within  the  experimental  error.  One  explanation  of  this 

asympoiotic  behavior  is  that  the  adsorption  sites,  which  could  be  a  small  fraction  of  the  total  surface 
sites,  e.g.,  kinks,  become  occupied  at  low  H2S  concentrations.  Thus,  the  effectiveness  of  H?S  in 

slowing  down  the  diffusion  of  H  atoms  on  the  surface  and/or  blocking  the  recombination  sites  with 
side  reaction  (6),  thereby  increasing  0h,  reaches  a  saturation  limit.  Another  explanation  is  that  H2S" 
adsorbs  readily  on  the  surface  and  at  rather  low  H2S  concentrations  in  the  solution  begins  to  interfere 
with  adsorption  sites  for  hydrogen  on  the  surface.  The  limit  of  this  effect  could  be  a  decrease  in  0h 
with  increasing  H2S  concentration,  which  is  not  observed  in  Figs.  8  and  9.  The  equilibrium  hydrogen 
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Figure  7  shows  that  the  plots  of  In  vs  i„  at  pH  =  2  are  linear  in  accord  with  Equation 
( 12)  for  C{},S=1C’^  M  and  ^,$-3x10*4  M.  But  for  the  lowest  HoS  concentration,  cJhS=2.5xlO'^  M, 

the  plot  nviv  not  be  completely  linear  and  is  consistent  with  the  h.e.r.  being  only  partially  activated,  as 
is  already  indicated  above  in  the  i_  vs  V ir  plot  (Figure  5(b))  which  is  linear  and  in  the  fj^  vs  i„  plot 

(Figure  6(b))  which  is  nonlinear.  Extended  analysis  showed  that  the  data  for  c„iS=2.5xl0‘^  M  are 

better  analyzed  with  the  F-T  corrected  I-P-Z  model. 

The  value  of  a  is  determined  from  Equations  (12)  and  (14)  by  an  iterative  procedure.  Using  an 
initial  value  of  a  =  0.5  and  the  k"  obtained  from  the  slope  of  the  (-y/ir  /i„)  vs.  i„  plot  (Eq.  (12)),  In 
(f;  :  )  is  known  for  various  values  of  ic  and  i„  and  as  a  function  of  T).  Thus,  a  regression  analysis 

of  Eq.  (14)  can  be  done  to  obtain  a  new  value  of  a.  Using  this  new  value  ot  a,  a  new  value  of  k"  is 

obtained  (from  Eq.  (12))  and  a  regression  analysis  is  again  performed  on  Eq.  (14).  This  procedure  is 
repeated  until  the  a  obtained  from  the  slope  of  the  In  (fjc  j  )  vs  T|  plot  (Eq.  (14))  converges  to  a  fixed 

value.  Then,  the  In  (fj  j^)  vs  tj  functions  are  plotted  (Figure  8),  and,  besides  a,  kg,  k"  and  i0'  are 

determined  from  the  regression  analysis  of  Equations  (12)  and  (14). 

Having  determined  k"  for  all  combinations  of  pH  and  cH<jS  (ie.  with  and  without  the  F-T 

correction  to  the  I-P-Z  model),  8pj  (calculated  using  Equation  (1 1))  vs  T|  can  be  plotted  and  these  are 
shown  in  Figures  9  (a)  and  (b).  These  olots  demonstrate  that  HoS  increases  the  hydrogen  surface 
coverage  (0pf)  quite  significantly  and  that  the  increase  is  asymptotic  with  increasing  c1JoS.  In  fact,  at 

pH  =  2  (Figure  9  (b)),  the  0j^  values  increase  significantly  with  small  additions  of  H^S,  but  do  not 
increase  any  further  with  higher  C|.  s  within  the  experimental  error.  One  explanation  of  this 

asympoiotic  behavior  in  that  the  adsorption  sites,  which  could  be  a  small  fraction  of  the  total  surface 
sites,  e.g.,  kinks,  become  occupied  at  low  H2S  concentrations.  Thus,  the  effectiveness  of  HgS  in 

slowing  down  the  diffusion  of  H  atoms  on  the  surface  and/or  blocking  the  recombination  sites  with 
side  reaction  (6),  thereby  increasing  0H.  reaches  a  saturation  limit.  Another  explanation  is  that  H2S' 
adsorbs  readily  on  the  surface  and  at  rather  low  H2S  concentrations  in  the  solution  begins  to  interfere 
with  adsorption  sites  for  hydrogen  on  the  surface.  The  limit  of  this  effect  could  be  a  decrease  in  0H 
with  increasing  H2S  concentration,  which  is  not  observed  in  Figs.  8  and  9.  The  equilibrium  hydrogen 
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coverage  (9e)  can  be  obtained  by  extrapolating  the  0j.j  vs  n  plots  to  rj  =  0,  and  these  plots  are  in 
Figure  9c.  Then,  i0  is  obtained  as:  i0=  \0*  (l-9e),  and  k°  is  obtained  from  Eq.  1. 

4.  HoS  Effects  on  the  Various  Kinetic  Parameters  and  the  h.e.r.  Transfer  Coefficient 

In  order  to  thoroughly  and  quantitatively  understand  the  mechanism  by  which  I-foS  enhances 
hydrogen  entry  into  iron,  alJ  of  the  important  quantities  are  plotted  versus  cH/(S  (Figures  10  to  15). 
Figure  10  shows  a  plot  of  i0  as  a  function  of  cHoS,  indicating  that  i0is  unaffected  by  I-^S.  On  the 
other  hand,  a  is  increased  significantly  by  small  additions  of  f-foS  as  shown  in  Figure  1 1.  Figure  1 1 
also  shows  within  the  experimental  scatter  that  a  asymptotically  reaches  a  concentration  independent 
value  at  higher  concentrations  of  HbS.  The  general  increase  in  a  can  be  explained  by  the  Kawashima, 
et.  al.,  bridge  mechanism  (described  earlier)  that  facilitates  easier  proton  discharge,  and  the  asymptotic 
behavior  of  a  with  cH>iS  may  be  a  consequence  of  the  H2S'  saturation  discussed  above. 

A  plot  of  the  discharge  rate  constant  (k  j°)  as  a  function  of  cH2S,  Figure  12,  shows  that  kj°  is 
unaffected  by  H2S.  The  large  scatter  in  Figure  12  may  be  due  to  the  actual  desorption  mechanism 
being  partly  chemical  recombination  and  partly  electrochemical,  especially  when  the  cH2$  is  high  since, 
as  is  reported  below,  chemical  recombination  less  readily  occurs  with  increasing  FbS  concentration. 
Due  to  the  extreme  complexity  of  this  problem  and  a  lack  of  the  pertinent  kinetic  data,  an  analysis 
considering  both  of  these  desorption  mechanisms  simultaneously  is  not  being  pursued  at  present. 

Another  important  parameter  is  the  thickness  dependent  absorption-adsorption  constant  (k") 
plotted  as  a  function  of  cH2S  (Figure  13),  showing  that  k"  is  not  altered  significantly  by  H2S.  The 
slight  increase  that  may  be  indicated  in  Figure  13  means  that  absorption  is  favored  over  adsorption 
with  increasing  H2S  concentration.  Hence  along  with  the  increasing  a  and  decreasing  k3  (reported 
below),  an  increasing  k"  may  also  contribute  to  the  increase  in  steady  state  permeation  with  increasing 
H2S  concentration.  An  increase  in  k"  could  be  caused  by  any  of  the  factors  associated  with  the  energy 
of  adsorpdon  becoming  coverage  dependent.  The  data  in  Figure  13  include  those  for  pH  =  1  (without 
F-T-c)  and  also  those  for  pH  -  2  (with  F-T-c)  and  they  seem  to  merge  quite  well,  further  supporting 
the  proposed  mechanism. 
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From  the  foregoing  discussion,  the  HiS  bridging  during  hydrogen  discharge  (reactions  (7)  & 
(8))  and  the  activation  of  the  h.e.r.  by  HjS  (ie.  the  slowing  down  of  the  diffusion  of  hydrogen 

adatoms  and/or  the  blocking  of  the  hydrogen  recombination  sites)  arc  considered  to  be  the  two  major 
effects  of  HoS.  A  further  confirmation  of  the  importance  of  these  effects  of  HpS  would  be  a 
significant  decrease  in  the  recombination  rate  constant  with  increasing  C|}oS.  Indeed,  that  this  occurs, 

is  clear  from  Figure  14,  which  shows  k3  decreases  with  addition  of  H2S  to  the  solution  and  as  with 
the  other  parameters,  then  levels  off  with  further  additions  of  IrtoS. 

5.  Hydrogen  Evolution  Reaction  Mechanism 

Finally,  the  potential  range  in  which  the  coupled  discharge-recombination  mechanism  operates 
can  also  be  determined  from  the  I-P-Z  model  ^)  as  follows: 

rig1  =[ln(10k1/k3)]/(aa)  (15) 

ncu  “  [In  (0-1  k1/k3)]/(aa)  (16) 

where  T|  *  is  the  lower  limit  and  ti  u  is  die  upper  limit  of  the  potential  range.  The  calculated  values  of 
c  c 

r|c1  and  ticu  arc  plotted  in  Figure  15,  as  a  function  of  cHzS  and  cH+.  The  two  solid  slanted  lines  are 

drawn,  one  through  the  lower  points  and  one  through  the  upper  points,  to  indicate  the  lower  and  upper 
bounds  of  the  hydrogen  overvoltages  calculated  for  the  coupled  discharge-recombination  mechanism 
using  the  ki,  k3  and  a  values  obtained  by  application  of  the  model  to  the  data  for  the  different  FbS 
concentrations,  while  the  horizontal  lines  at  T|  =-0.35  and  Tj  =-0.45V  indicate  the  lower  and  upper 
bounds  of  the  experimentally  measured  hydrogen  overvoltages.  For  example,  in  the  absence  of  H3S 

(which  corresponds  to  the  abscissa  value  of  zero),  the  experimental  overvoltage  range  is  at  less 

negative  electrode  potentials  than  the  overvoltage  range  for  operation  of  the  coupled  discharge- 

recombination  mechanism.  So,  the  discharge  reaction  may  be  rate  controlling  in  this  case.  But  with 
increasing  cHzS,  the  overpotential  range  for  operation  of  the  coupled  discharge-recombination 

mechanism  shifts  to  lower  overpotentials  so  that  above  a  certain  combination  of  the  concentrations  of 
H3S  and  H+  (i.e.,  >  0.2),  this  range  overlaps  the  overpotential  range  of  the 

experiments  (shaded  region  in  Figure  15).  Thus,  with  increasing  H2S  concentration,  the  h.e.r. 
mechanism  shifts  from  proton  discharge  controlled  to  coupled  discharge-recombination  controlled. 
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At  higher  pH  values,  one  would  expect  that  the  effect  of  l-foS  in  slowing  the  Effusion  of 

hydrogen  adatoms  will  be  even  larger  because  of  the  higher  (I  hSl/lH4-)  ratio.  What  this  would  mean 
is  that  the  mechanism  of  I-b  formation  tie.  hydrogen  f  .  ion)  may  have  to  citange  into  an 

electrochemical  desorption  mechanism  at  higher  pH  values.  '11ns  prediction  is  fully  consistent  with  the 
findings  of  Murayama  et  al  ^  9)  ;n  that  the  electrochemical  combination  route  lie.  Ha(j  +  H+  +  e‘  — > 
H2)  dominates  over  the  chemical  combination  route  for  99.9C&  iron  in  acetate  solutions  of  pH  >  3  in 
the  presence  of  HyS. 

One  final  comment  is  to  be  made  with  regard  to  the  importance  of  the  Equation  (3b)  in 
enlightening  our  understanding  of  the  surface  and  near  surface  reaction  kinetics.  Rewriting  Equation 
(3b),  we  obtain 

Ck'T1  =  OWW  +  (DA'abslf1-)'1  (17) 

Equation  (17)  has  been  discussed  clscwhcrcU^),  ulc  essence  of  Equation  (17)  is  that  if  the  plot  of 
1/k"  vs  1/L  is  linear,  the  values  of  kajjS  and  ka£js  are  obtained  (knowing  Dj).  Titus,  from  the  complete 

polarization  and  permeation  data  for  different  membrane  thicknesses  (L),  one  can  systematically 
investigate  the  important  surface  properties  relating  to  the  hydrogen  absorption  and  adsorption  rate 
constants. 

CONCLUSIONS 

1.  Cathodic  hydrogen  charging  experiments  were  performed  on  Ferrovac  E-Iron  membranes  in 
perchloric  acid  solutions  of  pH  1  and  2  with  and  witltout  fyS.  The  steady-state,  hydrogen  permeation 
current  density  (ioO  was  greatly  enhanced  in  the  presence  of  HjS,  while  the  hydrogen  overvoltage  (t|) 

was  decreased. 

2.  The  recently  developed  I-P-Z  model  (13)  is  applicable  to  the  h.e.r, /hydrogen  permeation 
process  in  the  absence  of  H2S  and  for  the  pH  1  solutions  of  lower  I-hS  concentrations.  From  the 

analysis  the  various  rate  constants,  exchange  current  density,  transfer  coefficient  and  hydrogen 
surface  coverage  were  obtained. 
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3.  For  the  higher  pH  and  higher  I-bS  concentrations,  the  I-P-Z  model  was  modifted  by  utilizing 
the  Frumkin-Temkin  correction  for  the  proton  discharge  and  hydrogen  recombination  reactions.  This 

correction  was  necessary  due  to  the  energy  of  adsorption  becoming  coverage  dependent,  probably  a 
result  of  the  product  of  the  fast  side  reaction,  ILS  +  c“  — >  l  hS'^,  restricting  the  diffusion  of 

hydrogen  adatoms  and/or  blocking  the  hydrogen  recombination  sites  on  the  iron  surface. 

4.  The  hydrogen  surface  coverage  (0j.j)  and  the  transfer  coefficient  (a)  arc  enhanced  while  the 

recombination  rate  constant  (k3)  is  decreased  in  the  presence  of  ILS.  The  increase  in  a  explains  the 
observed  decrease  in  ovcrpotcntial  (tj)  on  the  basis  of  the  Kawashima  ct  al  (3)  bridging  effect  of  IbS" 
tic.,  the  overall  reaction:  H30+  — >  Ha(js  +  l-foS  +  HoO). 

5.  On  the  basis  of  the  bridging  effect  and  the  slowing  of  the  recombination  reaction,  all  of  the 
observations  and  calculated  parametric  changes  by  HoS  have  been  rationally  explained.  These 

parametric  changes  are  the  increase  in  the  transfer  coefficient  (a)  and  the  decrease  in  the  recombination 
reaction  rate  constant  (kg).  In  addition,  the  overvoltage  range  for  operation  of  the  coupled  discharge- 

recombination  mechanism  is  predicted,  by  application  of  the  I-P-Z  model  to  the  data,  to  be  at  lower 
overpotentials  in  the  presence  of  H2S  and  to  overlap  the  experimentally  measured  overpotentials  above 
a  certain  IbS  concentration,  i.e.,  the  h.e.r.  mechanism  shifts  from  discharge  controlled  to  coupled 
discharge-recombination  controlled  with  increasing  I  hS  concentration. 

6.  It  is  noteworthy  that  one  can  now  calculate  the  important  hydrogen  absorption  and  adsorption 
rate  constants,  by  applying  the  I-P-Z  model  to  the  electrochemical  permeation  data  obtained  as  a 
function  of  membrane  thickness. 
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LIST  OF  SYMBOLS 

a  a  constant,  F/RT,  (volts)'1 

aH+  activity  of  hydrogen  ions  in  the  electrolyte,  dimensionless 


A  amperes 

b  a  constant,  L/FDj,  mol  (A  cm)'  * 

C|J+  hydrogen  ion  concentration,  10'PK,  mol  liter1 
°INS  HoS  concentration  in  the  electrolyte,  mol  liter** 

cs  surface  hydrogen  concentration,  mol  cm~3 
Dj  bulk  hydrogen  diffusion  coefficient,  cm-  s' 1 
Ee«  equilibrium  potential  for  the  h.e.r.,  mV  vs  SHE 

F  Faraday  constant,  96500  C  g*eq'  1 
f  a  constant  =  y/RT,  dimensionless 

fjc  j^  a  variable,  igC^^'00^  tyo-bi^"),  A  cm'2 

fj«  a  variable,  iccaocn,  a  cm'2 
ie  charging  flux  or  current  density,  A  cm'— 
i„  steady  state  permeation  flux,  A  cm’2 
i0  exchange  current  densitv.  A  cm'- 

i0'  io/(l-0c),  A  cm'2 

ir  steady  state  evolution  flux,  A  cm'2 

k  j°  discharge  reaction  rate  constant,  mol  (cm2  s)' 1 

k  j  discharge  reaction  rate  coefficient  (=  i0'  /F),  mol  (cm2  s)' 1 
ko  recombination  reaction  rate  constant,  mol  (cm2  s)* 1 
kai)S  absorption  rate  constant,  mol  (cm2  s)*1 
k;i(«s  adsorption  rate  constant,  cm  s' 1 

k"  thickness  dependent  absorption-adsorption  constant,  mol  enr 

L  membrane  thickness,  cm 

M  mol  liter'  ‘ 

R  gas  constant,  8.314  J  (g-mol  k)'1 

T  temperature,  K 

Greek  Symbols 

a  transfer  coefficient,  dimensionless 

q  overvoltage,  V 

qc  cathodic  overvoltage,  V 

qc*  lower  limit  of  t|c,  V 

qc11  upper  limit  of  qc,  V 

0pj  surface  hydrogen  coverage,  dimensionless 
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( 

0C  equilibrium  surface  hydrogen  coverage,  dimensionless 

Y  gradient  of  the  apparent  standard  free  energy  of  adsorption  with  hydrogen  coverage,  J  g-moH 


REFERENCES 

1.  HnS  Corrosion  in  Oil  and  Gas  Production.  A  Compilation  of  Classic  Papers,  R.  D.  Tuttle  and  R. 
S.~Kane,  eds.,  NACE  (1981). 

2.  J.  F.  Newman  and  L.  L.  Shrcir,  Corrosion  Science,  Vol.  9,  p.  631  (1969). 

3.  A.  Kawashima,  K.  Hashimoto  and  S.  Shimodaini,  Corrosion,  Vol.  32,  p.  321  (1976). 

4.  P.  Sury,  Corrosion  Science,  Vol.  16,  p.  879  (1976). 

5.  M.  Hashimoto,  E.  Sato  and  T.  Murata,  Hvdrogcn  in  Metals  (Conf.  Proc.),  JIMS-2  (1980),  p. 
209 

6.  B.  J.  Bcrkowitz  and  H.  H.  Horowitz,  J.  EIcctrochem.  Soc.,  Vol.  129,  p.  468  (1,982). 

7.  J.  O'M.  Bockris,  International  Conference  or.  Stress  Corrosion  Cracking  and  Hvdrouen 
Emhrittlemcnt  of  Iron  Base  Allovs.  Fimiiny,  France,  1973,  R.  W.  Stachlect  al  eds.,  NACE-5 
(1977),  p.  2S6. 

8.  Z.  A.  Iofa  and  F.  L.  Kam,  Zashchita  Metallov,  Vol.  10,  p.  17  (1974);  Protection  of  Metals,  Vol. 
10,  p.  12  (1974). 

9.  B.  E.  Boucher,  Proggcdings.Qfjiig_l:Qi!aiiimgrniuional  Congress, on  Mciallic.ConBSian,  NACE 
(1969),  p.  550. 

10.  J.  O'M.  Bockris,  J.  McBrcen  and  L.  Nanis,  J.  EIcctrochem.  Soc.,  Vol.  1 12,  p.  1025  (1965). 

1 1.  P.  W.  Bolmer,  Corrosion,  Vol.  21,  p.  69  (1965). 

12.  R.  N.  Iyer,  H.  W.  Pickering  and  M.  Zamanzadch,  Scripta  Met.,  Vol.  22,  p.  91 1  (1988). 

13.  R.  N.  Iyer,  H.  W.  Pickering  and  M.  Zamanzadch,  J.  EIcctrochem.  Soc.,  136. 7463  (1989). 

14.  M.  A.  Dcvanathan  and  L.  Stachurski,  J.  EIcctrochem.  Soc.,  Vol.  1 1 1,  p.  619  (1964). 

15.  S.  S.  Chaucrjce,  B.  G.  Atcya  and  H.  W.  Pickering,  Met.  Trans.,  Vol.  9A,  p.  389  (1978). 

16.  Private  Discussions  with  Konrad  Weil,  Technische  Hochschule  Darmstadt,  Darmstadt,  Federal 
Republic  of  Germany. 

17.  E.  Hart  and  M.  Anbar,  The  Hvdrnted  Electron.  John  Wiley  &  Sons,  Inc.,  New  York  (1970),  p. 
113. 

18.  E.  Gileadi  and  B.  E.  Conway,  Modem  Aspects  of  Electrochemistry.  No.  3,  J.  O’M.  Bockris  and 
B.  E.  Conway,  eds.,  Butterworths,  Washington  (1964),  pp.  347-442. 

19.  H.  Murayama,  M.  Sakashita  and  N.  Sato,  Hydrogen  in  Metals  (Conf.  Proc.),  JIMS-2  (1980),  p. 
297. 


16 


FIGURE  CAPTIONS 

Figure  l:  Stability  plots  of  HoS,  MS"  and  S“"  as  a  function  of  pH.  ^ 

Figure  2:  Schematic  of  the  apparatus  for  the  hydrogen  permeation 

experiment. 

Figure  3:  Hydrogen  overvoltage,  n ,  vs.  hydrogen  charging  current  density  ic, 
for  different  H2S  concentrations  at  (a)  pH=l  and  (b)  pH=2. 

Figure  4:  Hydrogen  overvoltage,  T),  vs.  steady  state  hydrogen  permeation  current 
density,  ioo,  for  different  H9S  concentrations  at  (a)  pH=l  and  (b)  pl*I=2. 

Figure  5:  Steady  state  hydrogen  permeation  current  density,  i^,  vs.  square  root  of 
hydrogen  recombination  current  density,  Vir,  for  different  H2S 
concentrations  at  (a)  pH=l  and  (b)  pH=2. 

Figure  6:  Relationship  between  fj^  and  i,*,  for  different  I-I2S  concentrations  at 
(a)  pH=l  and  (b)  pH=2. 

Figure  7:  Relationship  between  'll/loo  and  i*,  for  different  H2S  concentrations  at  pH=2 

Figure  8:  Hydrogen  coverage-corrected  ic  vs  hydrogen  overvoltage,  n,  for 
different  PfoS  concentrations  atpH=2. 

Figure  9:  Hydrogen  coverage,  Gj-j,  vs  hydrogen  overvoltage,  r\,  for  different  PbS 
concentrations  at  (a)  pH=l  and  (b)  pH=2. 

Figure  10:  I-I2S  effect  on  the  h.e.r.  exchange  current  density,  i0. 

Figure  1 1:  I-I9S  effect  on  the  h.e.r.  transfer  coefficient,  a. 

Figure  12:  H9S  effect  on  the  discharge  rate  constant,  kf. 

Figure  13:  H2S  effect  on  the  absorption-adsorption  constant,  k". 

Figure  14:  H2S  effect  on  the  hydrogen  recombination  reaction  rate  constant,  k3. 

Figure  15:  H2S  effect  on  the  hydrogen  overvoltage  range  (between  sloped  lines)  of 
the  coupled  discharge-recombination  mechanism.  Points:  Calculated 
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(Eqs.  15  and  16)  for  different  ki  and  k3  values  obtained  by  application 
of  the  I-P-Z  model  to  the  permeation  data.  Cross  hatched  area: 
potential  range  of  experimentation  overlaps  potential  range  of  coupled 
discharge-recombination  mechanism. 


currenl 

source 


OIOA)  U,  'a6ei|OAJ3AO  ua6ojpfiH 


iiriTi 


Ill 


•IU 


0.3 


(I/IOUJUJ) 

IS2HJ 


[H2S] 

(mmol/l) 


0.15 


pH  i0  (A/sq.cm.) 

1  1.18*I(T7 

2  5.47*I0~8 


o 

CM 

10 

o 

T~ 

to 

o 

o 

o 

o 

o 

o 

o 

o 

>0V  ‘>o  Ut  3SB3J0UI 


Fig.  1 


10' 


(03S/UU0)  0l>| 


Fig.  12 


Fig.  14 


BASIC  DISTRIBUTION  .1ST 


Technical  and  Summary  Reports 

1988 

Oraanization 

Cooies 

Organization  Cooi 

Defense  Documentation  Center 

Naval  Air  Prop.  Test  Ctr. 

Camerson  Station 

Trenton,  NY  08628 

Alexandria,  VA  22314 

12 

ATTN:  Library  1 

Office  of  Naval  Research 

Naval  Contruction  Battallion 

Dept,  of  the  Navy 

Civil  Engineering  Laboratory 

800  N.  Quincy  Street 

• 

Port  Hueneme,  CA  93043 

Arlington,  VA  22217 

3 

ATTN:  Materials  Oiv.  1 

Attn:  Code  1131 

Naval  Research  Laboratory 

Naval  Electronics  Laboratory 

Washington,  DC  20375  ' 

San  Oiego,  CA  92152 

ATTN:  Codes  6000 

ATTN:  Electron  Materials 

6300 

2627 

Sciences  Division  1 

Naval  Air  Development  Center 

Naval  Missile  Center 

Code  606 

Materials  Consultant 

Warminister,  PA  18974 

Code  3312-1 

ATTN:  Dr.  J.  DeLuccia 

Point  Mugu,  CA  92041  1 

Commanding  Officer 

Corrmander 

Naval  Surface  Weapons  Center 

David  Taylor  Research 

White  Oak  Laboratory 

Center 

Silver  Spring,  MD  20910 

Bethesda,  MD  20084  1 

ATTN:  Library 

1 

Naval  Oceans  Systems  Center 

Naval  Underwater  System  Ctr. 

San  Diego,  CA  92132 

ATTN:  Library 

Newport,  RI  02840 

1 

ATTN:  Library  1 

Naval  Postgraduate  School 

Naval  Weapons  Center 

Monterey,  CA  93940 

China  Lake,  CA  93555 

ATTN:  Mechanical  Engineering 

ATTN:  Library  1 

Department 

1 

Naval  Air  Systems  Command 

NASA 

Washington,  DC  20360 

Attn:  Code  310A 

Lewis  Research  Center 

1 

21000  Brookpark  Road 

Code  53048 

1 

Cleveland,  OH  44135 

Code  931A 

1 

ATTN:  Library  1 

Naval  Sea  System  Command 

National  Bureau  of  Standards 

Washinoton,  DC  20362 

Gaithersburg,  MD  20899 

ATTN:  Code  05M 

1 

Attn:  Metallurgy  Division  1 

Code  05R 

1 

Ceramics  Division  1 

Fracture  &  Deformation 
Division  1 


Htva 1  Facilities  Engineering 

Command 

Alexandria,  VA  22331 
ATTN:  Code  03 


Scientific  Advisor 
Commandant  of  the  Marine  Corps 
Washinoton,  DC  20380 
ATTN:  Code  AX 


1 


1 


Army  Research  Office 
P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
ATTN:  Metallurgy  &  Ceramics 

Program  1 

Army  Materials  and  Mechanics 
Research  Center 
Watertown,  MA  02172 
ATTN:  Research  Programs  O.ffice  1 

Air  Force  Office  of  Scientific 
Research/NE 
8uilding  410 
Bolling  Air  Force  Base 
Washington,  DC  20332 
ATTN:  Electronics  &  Materials 

Science  Directorate  1 


NASA  Headquarters 

Washington,  OC  20546 

Attn:.  Code  RM  1 


Defense  Metals  4  Ceramics 
Information  Center 
Battelle  Memorial  Inst. 

505  King  Avenue 

Columbus,  OH  43201  1 

Metals  and  Ceramics  Div. 

Oak  Ridge  National  Laboratory 
P.0.  Box  X 

Oak  Ridge,  TN  37380  1 

Los  Alamos  Scientific  Lab. 

P.0.  Box  1663 

Los  Alamos,  NM  87544 

ATTN:  Report  Librarian  1 


Argonne  .National  Laboratory 
Metallurgy  Division 
P.0.  Box  229 

Lemont,  IL  60439  1 

Brookhaven  National  Laboratory 
Technical  Information  Oivision 
Upton,  Long  Island 
New  York  11973 
Attn:  Research  Library  1 


Lawrence  Radiation  Lab. 
Library 

Building  50,  Room  134 
Berkely,  CA  1 

David  Taylor  Research  Ctr 
Annapolis,  MD  21402-5067 
ATTN:  Code  281  1 

Code  2813  1 

Code  0115  1 


RE/1131/88/75 


4315  (036) 

Supplemental  Distribution  List  Feb  1988 


Prof.  I.H.  Bernstein 
Illinois  Institute  of  Technology 
I IT  Center 
Chicago,  Ill  60615 

Prof.  H.K.  Birnbaun 
Dept,  of  Metallurgy  &  Mining  Eng. 
University  of  Illinois 
Urbana,  Ill  61801 


Prof.  H.W.  Pickering 
Dept,  of  Materials  Science  and  Eng. 
The  Pennsylvania  State  University 
University  Park,  PA  16802  v 


Prof.  D.J.-  Duquette 
Dept,  of  Metallurgical  Eng. 
Rensselaer  Polytechnic  Inst. 
Troy,  NY  12181 


Prof.  J.P.  Hirth 

Dept,  of  Metallurgical  Eng. 

The  Ohio  State  University 
116  West  19th  Avenue 
Columbus,  OH  43210-1179 

Prof.  H.  Leidheiser,  Jr. 

Center  for  Coatings  and  Surface  Research 
Sinclair  Laboratory,  Bid.  No.  7 
Lehigh  University 
Bethlehem,  PA  18015 

Dr.  M.  Kendig 

Rockwell  International  Science  Center 

1049  Camino  Dos  Rios 

P.0.  Box  1035 

Thousand  Oaks,  CA  91360 

Prof.  R.  A.  Rapp 

Oept.  of  Metallurgical  Eng. 

The  Ohio  State  University 
116  West  19th  Avenue 
Columbus,  OH  43210-1179 


Profs.  G.H.  Meier  and  F.S. Pettit 
Dept,  of  Metallurgical  and 
Materials  Eng. 

University  of  Pittsburgh 
Pittsburgh,  PA  15261 

Dr*.  W.  C.  Moshier 

Martin  Marietta  Laboratories 

1450  South  Rollino  Rd. 

Baltimore,  MD  21227-3898 

Prof.  P.J.  Moran 

Dept,  of  Materials  Science  &  Eng. 
The  Johns  Hopkins  University 
Baltimore,  MO  21218 

Prof.  J.  Kruger 

Dept,  of  Materials  Science  &  Eng. 
The  Johns  Hopkins  University 
Baltimore,  MD  21218 

Prof.  R.P.  Wei 

Dept,  of  Mechanical  Engineering 
and  Mechanics 
Lehigh  University 
Bethlehem,  PA  18015 


Prof.  W.H.  Hartt 
Department  of  Ocean  Engineering 
Florida  Atlantic  University 
Boca  Raton,  Florida  33431 


Dr.  B.G.  Pound 
SRI  International 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 

Prof.  C.R.  Clayton 
Department  of  Materials  Science 
&  Engineering 

State  University  of  New  York 
Stony  Brook 

Long  Island,  New  York  11794 


FVfrof.  Son’s  0.  Cahan 
gept.  of  Chemistry 
Case  Western  Reserve  Univ. 
Cleveland,  Ohio  44106 


Or.  K.  Sadananda 
Code  6323 

Naval  Research  Laboratory 
Washington,  D.C.  20375 

Prof.  M.E.  Orazem 
Dept,  of  Chemical  Engineering 
University  of  Virginia 
Charlottesville,  VA  22901 

Or.  G.R.  Yoder 
Code  6384 

Naval  Research  Laboratory 
Washington,  D.C.  20375 


Or.  N.  S.  Bornsteln 

United  Technologies  Research  Center 

East  Hartford,  CT  06108 


Dr.  A.L.  Moran 
Code  2812 

David  Taylor  Research  Center 
Annapolis,  MO  21402-5067 


Dr.  B.E.  Wilde 

Oept.  of  Metallurgical  Engineering 
The  Ohio  State  University 
116  West  19th  Avenue 
Columbus,  OH  43210-1179 


Prof.  G.R.  St.  Pierre 
Dept,  of  Metallurgical  Eng. 

The  Ohio  State  University 
116  West  19th  Avenue 
Columbus,  OH  43210-1179 

Prof.  G.  Simkovich 
Dept,  of  Materials  Science  &  Eng. 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dr.  E.  McCaffertv 
Code  6322 

Naval  Research  Laboratory 
Washington,  D.  C.  20375 

Dr.  J.A.  Sprague 
Code  4672 

.  Naval  Research  Laboratory 
Washington,  D.C.  20375 

Dr.  C.M.  Gilmore 
The  George  Washington  University 
School  of  Engineering  A  Applied 
Science 

Washington,  D.C.  20052 

Dr.  F.B.  Mansfeld 

Dept,  of  Materials  Science 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089 

Dr.  Ulrich  Stimming 
Dept,  of  Chemical  Eng.  &  Applied 
Chemistry 

Columbia  University 
Ndw  York,  N.Y.  10027 

Prof.  J.  O'M.  Bockris 
Dept,  of  Chemistry 
Texas  A  A  M  University 
College  Station,  TX  77843 


